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ABSTRACT 

Most of our knowledge of extrasolar planets rests on precise radial-velocity measure- 
ments, either for direct detection or for confirmation of the planetary origin of pho- 
tometric transit signals. This has limited our exploration of the parameter space of 
exoplanet hosts to solar- and later-type, sharp-lined stars. Here we extend the realm of 
stars with known planetary companions to include hot, fast-rotating stars. Planet-like 
transits have previously been reported in the lightcurve obtained by the SuperWASP 
survey of the A5 star HD15082 (WASP-33; V = 8.3, vsini = 86 km s'^). Here we 
report further photometry and time-series spectroscopy through three separate tran- 
sits, which we use to confirm the existence of a gas giant planet with an orbital period 
of 1.22d in orbit around HD15082. From the photometry and the properties of the 
planet signal travelling through the spectral line profiles during the transit we directly 
derive the size of the planet, the inclination and obliquity of its orbital plane, and its 
retrograde orbital motion relative to the spin of the star. This kind of analysis opens 
the way to studying the formation of planets around a whole new class of young, 
early-type stars, hence under different physical conditions and generally in an earlier 
stage of formation than in sharp-lined late-type stars. The reflex orbital motion of 
the star caused by the transiting planet is small, yielding an upper mass limit of 4.1 
M Jupiter on the planet. We also find evidence of a third body of sub-stellar mass in the 
system, which may explain the unusual orbit of the transiting planet. In HD 15082, the 
stellar line profiles also show evidence of non-radial pulsations, clearly distinct from 
the planetary transit signal. This raises the intriguing possibility that tides raised by 
the close-in planet may excite or amplify the pulsations in such stars. 

Key words: planetary systems - stars: rotation - stars: oscillations - binaries: eclips- 
ing - techniques: spectroscopic 
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1 INTRODUCTION 

More than 400 extra-solar planets have been discovered or- 
biting solar- and later-type stars, but very little is known so 
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far about planets orbiting intermediate-mass main-sequence 
stars 



Recent models of planet formation (e.g. Ida & Lin 
2005 Kennedy &: Kenyon||2008 1 suggest that gas-giant for- 
mation may peak among the late B stars. Radial-velocity 
searches for planets orbiting stars in this mass range become 
feasible only once the star has evolved into a late-type giant 
( Hatzes et aL][2005} [Johnson et "al]p007l [Sato et al.||2008| ), 
by which time any close-orbiting planets are likely to have 
been engulfed. Transit signals have been reported for main- 
sequence A-F-type stars, suggesting the possible presence 



of close-in planets (Christian et al. 2006 Lister et al. 2007 



Street et al.||2007| [Clarkson et al.||2007| [Kane et al.pOOSf , 



but measuring the small reflex motion that would confirm 
their planetary origin is generally not possible for these line- 
poor, fast-rotating stars. The hottest star known to passes a 
trans iting planet until now is OGLE2-TR-L9 ( Snellen et al. 
20091, which at spectral type F3 and vsmi = 33 km s~ is 



still not far from solar-type. As a result, we know very little 
yet about the frequency of giant-planet formation around 
earlier-type stars, with their harder radiation fields, more 
distant snow-lines and greater disc masses. 

Here we use the established technique of Doppler imag- 
ing to confirm the presence of a transiting planet around 
the bright, fast-rotating star HD 15082 (V = 8.3; Sp. type 
A5m, II sin i = 86 km s^^). The travelling spectral signature 
of a planet transiting the disc of the star during the light 
minimum yields information on the properties of the planet 
and its orbit. This complements or supplements that derived 
from standard photometric analyses, notably its size, retro- 
grade orbit, and non-alignment of the stellar and orbital spin 



2 DISCOVERY OF TRANSITS AND 
FOLLOW-UP PHOTOMETRY 

HD 15082 (WASP-33) is an early-type star in which flat- 
bottomed, planet-like transits recurring every 1.22 days were 



discovered by Christian et al. (2006| in the WASP survey 



(PoUacco et al. 20061. However, it rotates too rapidly to 



permit straightforward confirmation of the planet from pre- 
cise radial-velocity observations, and so was selected for the 
more sophisticated analysis reported here. 

Our first step was to follow up the WASP discovery with 
dedicated photometric observations of HD 15082 to improve 
the definition of the transit light curve. A partial transit was 
observed in the R band on 2006 November 13 using the CCD 
camera of the 0.95-m James Gregory Telescope ( JGT) at the 
St Andrews University Observatory. One complete and one 
partial _R-band transit were observed with the 60-cm tele- 
scope and CCD camera of the University of Keele ( |Pollacco| 
et al.||2008|) on 2006 December 11 and 2007 March 20. Fi- 



nally, complete transits were observed in R and I bands with 
the 35-cm Schmidt-Cassegrain telescope and CCD camera 



at the University of London Observatory, Mill Hill ( Fossey, 
Waldmann, fc Kipping"2009') on the nights of 2007 Novem- 



ber 14 and 2008 September 20. In all three instruments the 
unvignetted field of view includes several bright reference 
stars, which were used to derive the transit light curves 
shown in Fig. [l] using differential aperture photometry. The 
total transit duration is 2.72 hours from first to last contact; 



b 
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Figure 1. Transit light curves of HD 15082. Top to bottom: JGT 

2006 Nov 13, Keele 2006 Dec 11, Keele 2007 March 20, Mill Hill 

2007 Nov 14 (all K-band); and Mill Hill 2008 Sep 20 (/-band). 



ingress and egress last 16 minutes, and the transit depth is 
0.015 mag in both bands. 



3 EXPLORATORY SPECTROSCOPY 

Preliminary radial- velocity (RV) measurements of HD 15082 
were obtained at the Thiiringer Landessternwarte Tauten- 
burg (TLS) in order to determine the mass of the planet from 
the host star's orbital refiex motion. Observations were made 
during the periods 2006 November 30 - December 12 and 
2007 February 04 - 2007 March 05, using the coudeechelle 
spectrograph of the 2-m Alfred Jensch telescope at a resolv- 
ing power of A/AA = 67, 000, covering the wavelength range 
from 4700 to 7400 A. 

WASP-33 is part of the Tautenburg search for extraso- 
lar planets (e.g. [Guenther et al.|2009 |. Radial- velocity mea- 
surements are obtained using an iodine cell. As described 
in detail by _Hatzes et al.. (2005j , we fit an initial global 



wavelength solution to all orders of the spectrum of a ThAr 
lamp, then use the lines from the iodine cell to determine 
the instrumental profile (IP) and the instrumental shift si- 
multaneously with the observations. An accuracy o f 1.2 to 
1.7 m s"'^ is achieved for bright slowly-rotating stars (Hatzes 
fc Zechmeister|2007| ). 

In the case of WASP-33, however, the reduced stellar 
spectra show extreme rotational broadening, with usini ~ 
90 ± 10 km s~^. This precludes the direct determination 
of precise radial velocities relative to the iodine spectrum. 
Instead, relative radial velocities were computed by cross- 
correlation with the template spectrum of the star outside 
the iodine spectral region. 

The results are listed together with the signal-to-noise 
ratios (SNR) of the individual spectra in Table [l] and the 
radial-velocity curve is plotted in Fig. [2] There is a long- 
term trend in the radial velocities over the 95-day span of 
these observations, corresponding to a radial acceleration of 
— 18 ± 3 m s~^ d~^ in the centre of mass of the transiting 
pair. Although most of the stars in the Tautenburg radial- 
velocity surveys do not exhibit trends at this level ( |Guenther| 
et al.|2007[ ), we are cautious about attributing this trend to 
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Table 1. Radial velocities measured with the 2.0-m telescope at 
TLS. 
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another body in the system, because of the star's non-radial 
pulsations. 

The precision of the RV measurements is insufficient 
for a significant detection of the star's refiex motion on 
the 1.22-day orbit of the transiting b component from these 
data alone. The RMS scatter about the best-fitting solution 
is 0.99 km s~^, whereas the mean formally-estimated un- 
certainty on an individual RV measurement is 0.8 km s^^. 
The uncertainty in the measured radial-velocity amplitude 
yields a 99.9 percent upper limit of 4.1 Mjup for the mass 
of companion b, so if the transits are genuinely caused by a 
body orbiting the host star, it has to be of planetary rather 
than stellar mass. However, the possibility remained that 
the transits could be caused by a faint, spatially-unresolved 
eclipsing-binary companion. 



4 STELLAR PROPERTIES OF HD 15082 



[Grenier et aLHT999l l classify HD 15082 as an A5mA8F4 star, 
i.e. A5 from the Ca ll K line, A8 from the H lines, and F4 
from the metal lines. This would suggest that HD 15082 is a 
classical Am star with overabundances of iron-group metals 



and underabundances of Ca and Sc (Wolff 19831. Stromgren- 
Crawford uvbyfi photometry (Hauck & Mermilliod 19971 
yields T^s = 7430 K, logg = 4.21 and [M/H] = 0.21 via 
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Figure 2. Radial-velocity curve of HD 15082 obtained with the 
iodine cell at Tautenburg. The upper panel shows the long-term 
radial acceleration. The lower panel shows the phased radial ve- 
locities with a linear trend of 0.018 km s^^ d~^ subtracted. The 
dashed curve represents the best-fitting circular-orbit solution. 
The zero-point of the radial-velocity scale is arbitrary. 

from the General Catalogue of Photometric Datsj^ yields 
T^ff = 7471 ± 63 K, logg = 4.35 ± .07 and [M/H] = 



0.08 ± 0.09 with the KunzU et al. ( 19971 calibration. 



The template spectrum (taken without the iodine cell) 
was analysed using the uclsyn spectral synthesis pack- 
age (Smith||1992| ISmalley, Smith, & Dworetsky||2001l). A 



reasonable fit to the spectrum is obtained with Te// — 
7400 ± 200 K, logg = 4.3 ± 0.2 and [M/H] = 0.1 ± 0.2, in 
good agreement with the values derived from photometry. 
No obvious Am characteristics are visible in this spectrum 
other than slightly weak Call HSzK lines; a definitive anal- 
ysis of this issue will require a spectrum with much higher 
signal-to-noise. The stellar properties are summarized in Ta- 
bled 



5 TIME-RESOLVED TRANSIT SPECTRA 

During a transit, a planet blocks a small part of the stel- 
lar photosphere. The missing starlight no longer contributes 



the calibration of Smalley ( 1993 1, while Geneva photometry 



^ GCPD, http://obswww.unige.ch/gcpd/gcpd.html 
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Table 2. Stellar parameters of WASP-33 from Tautenburg tem- 
plate spectrum and catalogue data. 



Parameter 


Value 


RA (2000.0) 


02m26m51.06s 


Dec (2000.0) 


+37°33'01.7" 


V mag. 


8.3 


Spectral Type 


kA5hA8mF4 


Distance 


116 ± 16 pc 




7430 ± 100 K 


logs 


4.3 ± 0.2 


V sin i 


90 ± 10 km s-i 


[M/H] 


0.1 ± 0.2 



Notes: Spectral type from 
Distance from Hipparcos 



Grenier et al.|||1999|l 
( Ferryman et al.|1997l l 



to the overall spectral line profile, but appears as a re- 
solved "bump" within the rotational broadening function, 
in the same way as dark star spots distort the line profiles 
of a rotating star (Vogt & Penrod 1983^ If such a signa- 



ture is present and matches the transit geometry deduced 
from photometry, the presence of a planet may be confirmed 
even if the host star's reflex orbit cannot be reliably mea- 



sured. Jenkins et al. (2010) have also noted recently that the 
Rossiter effect provides a powerful method for distinguish- 
ing planets from blends in cases where classical methods of 
confirmation cannot be used. 

The spectral signature of a blended stellar binary would 
look very different. To produce a transit with a fractional 
depth of 0.01, even an object with a very deep eclipse should 
contribute at least one percent of the total system light. At 
a period of 1.2 days, the radial-velocity amplitude of the 
primary star will be of order 100 km s^^, giving a velocity 
drift of order tens of km/sec within the 3-hour duration of 
a typical transit. When the light of the brighter foreground 
star is subtracted, such a background signature is easily de- 
tectable in the residual line profiles. The Rossiter distortion 
will be imposed on the rotation profile of this drifting stellar 
signature, not on that of the foreground star. 

A new set of echelle spectra was therefore obtained at 
TLS on the night of 2008 October 18. Sixteen spectra with 
600s integration time were obtained at _R = 40, 000 from 
4708 to 7004 Awithout the iodine cell, starting just before 
ingress and ending after egress. A second dataset was ob- 
tained with the CS23 echelle spectrograph on the 2.7-m 
Harlan J. Smith Telescope at McDonald Observatory during 
the night of 2008 November 11/12, with a spectral resolv- 
ing power of i? = 60, 000 from 3837 to 10209 A and expo- 
sure times of 900s. This series also covers the entire transit. 
A third sequence with 540-s time resolution was obtained 
with the FIES spectrograph at the Nordic Optical Telescope 
(NOT) on La Palma on 2009 December 8, extending for 
~ 2 hours before and after the transit to characterise better 
the stellar pulsations seen earlier and discussed below. The 
wavelength range was 3911 to 7271 A with R = 67000. 

For each dataset, we combined the information in all the 
spectral lines in each echellogram into a single, composite 
profile with a high signal-to-noise ratio using Least-Squares 



Deconvolution (LSD). This method, pioneered by Donati et 
|al.| ( |1997[ ), computes the line profile which, when convolved 
with a forest of appropriately-weighted delta functions at 
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Figure 3. Least-squares deconvolved profiles of HD 15082 during 
the planetary transit of 2009 December 08 observed at the NOT. 
Data are shown as points with error bars. The limb-darkened 
stellar rotation profile is shown as a solid line. Time increases 
vertically. The planet signature is the positive bump that migrates 
from ~ 5 km s~^ in the first (bottom) profile to d ~ —40 km s~^ 
in the last (top) profile shown. 



the wavelengths of known spectral lines, gives an optimal 
inverse variance-weighted fit to the data. The velocity scale 
is tied to the reference frame of the Earth's atmosphere along 
the observer's line of sight using a separate LSD calculation 
based on a mask of narrow telluric absorption lines. This 
compensates for any drift in the instrumental profile to an 
accuracy of order 10 to 20 m s"^ This is more than sufficient 
for tracking features in the line profile that move through 
several tens of km s~^ during the transit. The signal-to-noise 
ratios of the composite profiles are 1070 for the TLS data, 
1730 for the McDonald data, and 1390 for the NOT data. 
The resulting time-series spectra during the NOT transit are 
shown in one-dimensional form Figure |3] The trailed spec- 
tra for all three data sets are shown in grey-scale form in 
Figure [4j with the mean line profile of the star itself sub- 
tracted. They exhibit a clear Doppler signature of a planet 
moving in the retrograde direction across the stellar rotation 
profile. They also show the characteristic pattern of stellar 
non-radial pulsations carried in the prograde direction by 
the stellar rotation. We discuss the planetary signal further 
in the following and the pulsations in Section [7| below. 

Fig. [4] shows the wealth of direct information on the 
planet's motion that is available from these data: First, 
ingress (bottom right) clearly occurs over the receding limb 
of the star, egress (top left) over the approaching limb - 
immediate evidence that the orbital motion of the planet is 
opposite to the rotation of the star, i.e. retrograde. The ve- 
locity width of the planet signal relative to that of the stellar 
line profile corresponds to the length of the chord traversed 
by the planet relative to the diameter of the stellar disc, 
i.e. the impact parameter of the transit, related to the or- 
bital inclination. Similarly, ingress clearly occurs closer to 
the spin axis of the star than the egress, proof that the orbit 
is inclined relative to the stellar equator as well as to the 
line of sight, and that the trace is not due to a spot on the 
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Figure 4. Time series of the residual average spectral line profile of HD 15082 during the transits of 2008 October 18 at Tautenburg 
(upper left); 2008 November 11/12 at McDonald (upper right); and 2009 December 8 at the NOT (lower left). Wavelength or radial 
velocity increases from left to right, time from bottom to top. The limb-darkened rotation profile of the star itself has been subtracted 
between the dashed lines at ±vsini, and the velocity of the star is denoted by a dotted vertical line with the four contact times of the 
transit indicated by " +" symbols. The planetary signal is the bright feature that moves from right-to-left coincident with the duration 
of the photometric transit. In the lower-right panel, the model planet signature has been subtracted from the NOT data, leaving only 
the pattern of non-radial pulsations. 



surface of the star (the transit speed confirms this as well). 
Finally, the width and depth of the planet trace are direct 
indicators of the radius and surface area of the planet, when 
allowance is made for the effects of the PSF, natural line 
width, and acceleration of the planet during each exposure. 



Information on the rotation of the star and the direc- 
tion of its spin axis relative to that of the orbit can also 
be derived from accurate radial-velocity data through the 
Rossiter-McLaughlin effect (e.g. [Anderson et al.|20ld| , but is 
available in a far more direct and less model-dependent way 
from the line-profile analysis discussed here, when the stel- 
lar rotation is large enough. The radius of the star and the 
inclination of the orbit can be determined both by the stan- 
dard analysis of the transit light curve and by the Doppler 
imaging performed here, which allows at least a consistency 
check. Which of the methods yields the more accurate re- 



sults depends on the system in question and the quality of 
the two datasets. 



6 COMBINED ANALYSIS 

In this section we integrate all the available data into a com- 
bined analysis to determine the properties of the planet and 
the star. This includes the WASP data from 2004 discussed 
by ( Christian et al.||2006 l plus a further 3658 observations 
obtained by the WASP survey in 2007. 

For the light curves, we modelled the transit geome- 
try and orbital parameters with the Markov-chain Monte- 
Carlo parameter fitting code described in detail by [CoUierl 
Cameron et all (|2007|) and iPoUacco et al.l (120081). Five pa- 



rameters are determined: the epoch To of mid-transit, the 
orbital period P, the planet/star radius ratio Rp/Rs, the 
scaled stellar radius Rs/a (where a is the orbital semi-major 
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Table 3. System parameters of HD 15082 as derived from our photometry, radial velocities and transit spectroscopy. 



Parameter 



Phot+RV 



TLS 



McD 



NOT 



Epoch of mid-transit, Tc (HJD) 
Orbital period, P (days) 
Planet/star radius ratio, Rp/Rs 
Scaled stellar radius, Rs/a 
Impact parameter b 
Centre-of-mass velocity, 7 (km s""*^) 
Radial acceleration, d'y/dt (km s~^ d^^) 
Radial- velocity amplitude, Ks (km s~^) 
Stellar rotation speed, vsini (km s~^) 
Projected obliquity, A (degrees) 
Intrinsic linewidth, "Ufwhm (km s~^) 
Stellar density, ps/pQ 
Stellar mass, M^/Mq 
Stellar radius, Rs/Rq 
Planet mass, Mp/AIj^p 
Planet radius, Rp/Rjup 
Orbital separation, a (AU) 
Orbital inclination, i (degrees) 



2454163.22373 ± 0.00026 
1.2198669 ±0.0000012 
0.1066 ±0.0009 
0.2640 ± 0.0057 



0.155 



+0.100 
-0.120 



-0.0183 ±0.0035 
< 0.59 



0.497 ± 0.024 
1.495 ± 0.031 
1.444 ± 0.034 

< 4.1 
1.497 ±0.045 
0.02555 ±0.00017 
87.67 ±1.81 



0.1072 ±0.0006 
0.2673 ±0.0015 
0.218 ±0.008 
-2.19 ±0.09 



86.04 ±0.13 
251.2 ± 1.0 
16.2 ±0.5 
0.472 ± 0.008 
1.497 ± 0.029 
1.469 ±0.013 

1.533 ±0.016 
0.02556 ±0.00017 
86.67 ± 0.13 



0.1071 ±0.0006 
0.2653 ±0.0016 
0.176 ±0.010 
-3.69 ±0.09 



85.64 ± 0.13 
254.2 ± 1.2 
19.2 ±0.6 
0.483 ± 0.008 
1.495 ±0.029 
1.458 ±0.013 

1.519 ±0.016 
0.02555 ±0.00017 
87.33 ±0.15 



0.1073 ±0.0006 
0.2661 ± 0.0015 
0.203 ±0.007 
-2.11 ± 0.05 



86.48 ±0.06 
251.6 ±0.7 
18.1 ± 0.3 
0.479 ± 0.007 
1.500 ± 0.030 
1.464 ± 0.013 

1.527 ± 0.015 
0.02557 ± 0.00014 
86.91 ±0.11 



axis) and the impact parameter b = a cos i/Rs, where i is 
the orbital inclination. Model hght curves are computed us- 



ing the analytic formulation of Mandel & Agol ( 2002 1 in 



the small-planet approximation with non-linear limb dark- 
ening. Slow small-scale photometric trends are removed in 
the process. 

The stellar density p is derived from Rs/a and P via 
Kepler's third Law ( [Seager fc Mallen-Ornelas|2003[ ): 



Ps 



47ri?i GP^ 



(1) 



The stellar mass is estimated at each step from the empir- 



ical calibration of Torres, Andersen, & Gimenez ( 2009 1 as 
modified by |Enoch et al.| ( |2010[ ). The radius of the star, and 
hence that of the planet, then follow. 

A circular orbit with a constant radial-acceleration term 
was fit to the radial-velocity data to derive the stellar reflex 
velocity amplitude Kg. At each step in the Markov chain, 
the centre-of-mass velocity 7 at the fiducial epoch of transit 
is computed as the mean of the residuals. 

The MCMC solution combining the photometry and 
the Tautenburg radial velocities is given in the first column 
of Table |3] The secular trend in the values of 7 from 2006 
and 2007 is found to be -18 ± 3 m s"^ d"\ If attributed 
to a binary companion (denoted as component c), the de- 
rived mass would be of order Mc = 0.036{a^/AUf Mq. 
The outer period appears to be at least 200 days, suggest- 
ing that component c, if real, should have a mass at least 
20 times that of Jupiter. This interpretation should, how- 
ever, be treated with caution. The non-radial pulsations dis- 
tort the line-profile shapes strongly, introducing substantial 
radial- velocity shifts which could introduce a spurious radial 
acceleration into the orbital solution. The introduction of a 
linear trend into the model reduces the RMS scatter about 
the fitted solution from 1.30 km s^^ to 0.99 km s^^. An 
F-test on the residuals about the fitted solution with and 
without the linear trend returns a 17 percent probability 
that their variances are not significantly different. 



For the line profile analysis we incorporated each of the 
three time-resolved spectral sequences in the MCMC solu- 
tion. We used the profile decomposition method of [Collierl 



Cameron et al. ( 2010 1 to determine the rotation speed of the 



star and the projected orientation of the orbital and stellar 
spin axes. The model line profile is a limb-darkened rota- 
tion profile convolved with a Gaussian corresponding to the 
intrinsic photospheric line profile plus instrumental broad- 
ening. The planet's Doppler shadow is a travelling Gaus- 
sian of the same width, computed as described by |Man-| 
del & Agol (20021. The changing form of the line profile 



is thus modelled with five additional MCMC parameters 



(Collier Cameron et al. 20101: the stellar rotation speed 



wsini, the difference A between the position angles of the 
orbital and stellar rotation axes in the plane of the sky; 
the impact parameter 6 of the planet's trajectory; the width 
'i'FWHM of the Gaussian intrinsic profile; and a nightly zero- 
point correction 7 to the centre-of-mass velocity at the time 
of observation. The velocity difference between the planet 
shadow at ingress and egress relative to the stellar veloc- 
ity provides a better constraint on the value of the impact 
parameter b than the photometry alone, as noted above: 
Avi^c ~ bsin A ± V 1 — fc'^ cos A. 

Individual solutions are given for each set of transit 
spectroscopy in Table|3] combined with the photometry and 
Tautenburg radial velocities. The system centre-of-mass ve- 
locity differs slightly when different sets of transit spectra 
are used. This is partly attributable to the outer orbit, and 
possibly to zero-point differences between the different in- 
struments. The non-radial pulsations, however, are likely 
to be the dominant additional source of error in the pro- 
file fits. They introduce line-profile distortions which can- 
not be accounted for in the modelling procedure, and whose 
amplitudes approach that of the planet's Doppler shadow. 
The parameters b, wsini, A and ufwhm described in the 
paragraph above are the ones most directly affected by the 
pulsations, and indeed their fitted values disagree by sev- 
eral sigma among the three individual sets of observations. 
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Figure 5. The effective temperature and inverse cube root of 
density of HD 15082 are shown in relation to the evolutionary 
tracks of [Girardi et al.] l |2000[ | for stars of mass 1.5, 1.6, 1.7 and 
1.8 Mq, at ages of 63 Myr (solid), 320 Myr (dot-dash), 630 Myr 
(triple dot-dash), 795 Myr (dashed) and 1 Gyr (dotted). The error 
bar on the stellar density is smaller than the plot symbol. A mass 
of 1.55 ±0.04 M0 and an age less than 25 Myr are indicated. The 
tracks and isochrones are computed for [Fe/H]=0.08, correspond- 
ing to heavy-element abundance Z=0.0228. Decreasing [Fe/H] by 
0.09 dex lowers the zero-age main-sequence value of {p/ Pq)~^^^ 
by 0.03 for the youngest isochrone, and raises T^ff by 150K for a 
1.5Mq stellar model. 



Figure 6. The effective temperature and inverse cube root of 
density of HD 15082 are shown in relation to the Yale-Yonsei 
isochrones ( [Demarque et al.|[2004[ | for stars of [M/H] = -0.10 
(dot-dashed lines), 0.00 (dotted lines) and 0.10 (solid lines). The 
isochrone ages are 50, 100 and 400 Myr (bottom to top) in each 
case. The observed stellar density is clearly inconsistent with the 
upper part of the error range on [M/H] quoted in Table [5] the 
tracks shown here match the range from the nominal [M/H]=0.1 
to the one-sigma lower limit [M/H]=-0.1. If the star has solar 
metallicity, it is about 90 Myr old. If it it 0.1 dex sub-solar, the 
age increases to 300 Myr. Given the lower limit on the metallicity 
derived from the stellar spectra, we infer an upper age limit of 
~ 400 Myr. 



They propagate directly into the derived value for the or- 
bital inclination, which also exhibits a scatter greater than 
expected from the formal errors retmned by MCMC. For 
these parameters, the scatter among the different columns 
is a better indicator of the true uncertainty than the formal 
errors associated with the individual values. For all other 
quantities, the mutual agreement is excellent. 

Even including the uncertainties introduced by the pul- 
sations, the impact parameter is very well determined by the 
trailed spectra. With the help of the photometry, the stellar 
density is then determined to a precision of order 2 percent. 
This allows the stellar age and mass to be estimated in re- 
lation to isochrones and evolutionary tracks p lotted in the 
(Toff,(p./p0)"^/^) plane ( [Sozzetti et al.|2007| ). In Fig.[5]we 
show the location of HD 15082 in relation to the evolutionary 
models of [Girardi eTah] pOOO| for stars with [Fe/H]=0.08 
and masses in the range 1.5 to 1.8 Mq. The stellar mass is 
found to be 1.55 ± 0.04 Mq, in good agreement with the 
value 1.50 ± 0.03 Mq obtained from the empirical calibra- 
tion built into the MCMC analysis, as listed in Table |3] If 
we use evolutionary tracks interpolated for [M/H] =-0.1, we 
obtain a stellar mass of 1.50 Mq, which is slightly closer to 
the value obtained from the eclipsing-binary calibration. 

The stellar density is seen to be very close to the main- 
sequence value for stars in this mass range. Comparison with 



the Yonsei-Yale isochrones of Demarque et al. ( 2004 1 shows 



HD 15082 to be in the early part of its main-sequence life, 
with an isochrone age no greater than 400 Myr or so after 
making allowance for a 0.2 dex uncertainty in the metallicity 
(Fig.[6|. It is too dense to lie on or above the main sequence 
for metallicities greater than solar. 

Assuming a stellar TcS = 7430 ± lOOK, zero albedo and 



uniform redistribution of incident flux, the equilibrium tem- 
perature of the planet is 2710 ± 50 K. 



7 NON-RADIAL PULSATIONS IN THE HOST 
STAR 

Removing the planet's Doppler shadow from the trailed 
spectra highlights the underlying pattern of non-radial pul- 
sations of the star. Note that the radial component of the 
pulsations vanishes near the centre of the line profile (Fig. |4] 
lower right panel), a behaviour indicative of g- modes. Over- 
all, the pattern appears to be that of a sectoral non-radial 
mode with £ = 4 ± 2, suggesting that the star is a 7 Dor- 
type variable. The pulsation amplitude appears greater in 
the blue wing of the profile than the red. This phenomenon 
was noted earlier by e.g. Uytterhoeven et al. (2008) who 



studied HD 49434, a star remarkably similar in Tcs and 
wsini to HD 15082. HD 49434 exhibits several pulsation 
modes with frequencies similar to the ~ 4 cycles/day tidal 
forcing frequency imposed on HD 15082 by its planet. It has 
been suggested ( |Willems fc Aerts 2002 1 that stellar pulsa- 
tion modes may be excited by resonant tidal forcing by a 
close binary companion, and [Fekel, Warner, fc Kaye| ( |2003[ ) 
reported that a surprisingly high fraction of 7 Dor pulsators 
have close binary companions. Similar resonances are in- 
voked in the study of dissipative tidal spin-orbit coupling in 
hot- Jupiter systems ( Ogilvie fc Lin|[2007 Barker & Ogilvie 
2009 1 . Whether a retrograde planet can excite pulsations in 



a star rotating in the opposite direction remains an open 
question, however. 
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8 SUMMARY AND CONCLUSIONS 

Our time-series spectroscopy has provided the first confir- 
mation of a gas-giant planet transiting a rapidly-rotating 
A-type main-sequence star, proving that A stars do occa- 
sionally harbour close-in planets. The result opens up a new 
part of the parameter space of exoplanet host stars, where 
planets are formed under different physical conditions than 
around solar-type stars, and on shorter time-scales, ft will 
be interesting to compare the properties of planets orbit- 
ing young and intermediate-age stars as determined through 
the Doppler imaging technique to study the early dynamical 
evolution of planetary orbits in different environments. 

With a radius of 1.46 Rjup HD 15082b is one of the more 
bloated exoplanets yet discovered. The large radii of some of 
the hot Jupi ters, notably TrE S-4b ( [Mandushev et al. 12007]), 
W ASP- 1 2b ( |Hebb et al.|2009l ) and WASP- 17b ( [Anderson et 
|al.||2010| ) remain a challenge to theories of planetary struc- 
ture. Some of the leading contenders are dissipation of tidal 
energy in the interior of a planet during orbit circularisa- 
tion ( [Bodenheimer, Lin, fc Mardling||2001[ ), advection of a 
small fraction of incident stellar flux into the deep interior 
I Guillot & Showman 2002 1, the presence or absence of a sub- 
stantial rock/ice core ( [Fortney, Marley, fc Barnes|[2007j ) or 
some combination of the above, with a possible dependence 
of core mass on host-star metallicity ( Fressin et al.|2007[ |. 

Without a measurable radial- velocity orbit, the mass 
and density of the planet remain unknown, making it diffi- 
cult to estimate either its core mass or its orbital eccentric- 
ity. In its 1.2-day orbit, however, HD 15082b is the most 
strongly-irradiated planet yet found, with an equilibrium 
temperature of order 2710 K assuming isotropic re-radiation 
and a low albedo. This may be at least partly responsible for 
the planet's large radius (e.g. Guillot & Showman (2002|). 



At V = 8.3, HD 15082 is also the fifth-brightest star yet 
found to host a transiting planet. This will make it an at- 
tractive target for future infrared secondary-eclipse studies 
IjCharbonneau et al.||2005{ [Deming et al.||2005| [Knutson et| 
al. 20071. In particular, the displacement of the secondary 



eclipse from phase 0.5 and the duration of the secondary 
eclipse will allow the orbital eccentricity and argument of 
pericentre to be determined ( Carter et al.|2008 l. 

HD 15082b joins the growing ranks of retrograde plan- 



ets, notablyWAS]M7b^^ and HAT-P- 

9|). These and other 



7b ( |Winn et alT]|2009| [Narita et al.||20"0^ 

strongly misaligned planets present a challenge to theories of 
disc migration (e.g. Lin, Bodenheimer, fc Richardson|1996 1. 
Since the star and the planets form from the same disc, 
their respective spin and orbital angular momentum vec- 
tors should be aligned. The misaligned planets are, however, 
plausibly explained by alternative theories involving com- 



binations of planet-planet scattering (Rasio & Ford 1996 
Weidenschilling & Marzari 1996J, Kozai oscillations (Kozai 



1962 



Wu fc Murray 2003) and tidal damping (Eggleton & 



Kiseieva-Eggleton||2001. .Fabrycky fc Tremaine||2007| |Naga- 



sawa, Ida, fc Bessho 2008 1. A key future test of these theories 



will involve searching for additional, exterior bodies in these 
systems, which might have excited Kozai oscillations in the 
past. Continued radial- velocity monitoring of this star, and 
a deep search for faint common-proper-motion companions, 
would therefore be of interest. 

The host star is a non-radial pulsator, apparently of 



the 7 Dor type. A full study of the pulsation modes in HD 
15082 is beyond the scope of the present paper, but could 
in the future provide independent confirmation of the stellar 
density derived here from the transit geometry. 
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